Quasiperiodic events in an earthquake model 



0. Ramos 1 , E. Altshuler 2 , and K. J. Makty 1 
1 Department of Physics, University of Oslo, P.O.B. 1048, Blindern N-0316, Oslo, Norway 
"Henri Poincare" Group of Complex Systems, Physics Faculty, University of Havana, 10400 Havana, Cuba 

(Dated: February 2, 2008) 

We introduce a modification of the OFC earthquake model [Phys. Rev. Lett. 68, 1244 (1992)] in 
order to improve resemblance with the Burridge and Knopoff mechanical model and with possible 
laboratory experiments. A constant force continually drives the system, and thresholds are dis- 
tributed randomly following a narrow distribution. We find quasiperiodic behavior in the avalanche 
time series with a period proportional to the degree of dissipation of the system. Periodicity is not 
as robust as criticality when the threshold force distribution widens; and foreshocks and aftershocks 
are connected to the observed periodicity. 
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The concept of Self-organized criticality (SOC) Q, in- 
troduced by Bak, Tang and Wiesenfeld in 1987, was an 
attempt to explain the appearance of scale invariance in 
nature. In their sandpile model both the random, slow 
addition of "blocks" on a two dimensional lattice and a 
simple, local, and conservative rule drive the system into 
a critical state where power law distributed avalanches 
maintain a steady regime far from equilibrium. There 
is no correlation between the avalanches, and eventu- 
ally they reach the boundaries of the lattice liberating 
the excess of energy. Five years later Olami, Feder and 
Christensen (OFC) made an important contribution to 
the SOC ideas by mapping the Burridge and Knopoff 
spring-block model into a nonconservative cellular au- 
tomata simulating the earthquake's behavior and 
introducing dissipation in the family of SOC systems. 
The fact that avalanches are uncorrelated in the sandpile 
model has been used as an argument to propose that it 
is not possible to predict real earthquakes |g. However, 
foreshocks, aftershocks and clustering properties @ indi- 
cate the existence of correlation between different events. 
Many seismologists believe that large earthquakes are 
quasiperiodic [7J, |8| , but periodic behavior has appeared 
in theoretical models only as a special or as a trivial so- 
lution [Hnnm. 

The spring-block model consists in a two dimensional 
array of blocks on a flat surface. Each block is connected 
with its four nearest neighbors, and in the vertical di- 
rection, to a driving plate which moves horizontally at 
velocity v. The connections are made by springs, and 
when the force acting on a block overcomes the static 
friction with the surface, the block slips. Then a redistri- 
bution of forces takes place in the neighbors that even- 
tually trigger new displacements. In the OFC model the 
force on a block is stored in a site of a squared lattice, 
and the static friction threshold has the same value for 
all blocks. Starting from a random distribution of forces, 
the site closest to the threshold is found and the exact 
force necessary to provoke a slip in this block is added 
to every site of the grid. This infinitely accurate tuning 
is only possible in the mechanical model if the displace- 
ment of the plate is infinitely slow (v — > 0), considering 



the fact that the real time resolution is finite. When a 
site reaches the threshold it is set to zero, and a frac- 
tion a of its force is redistributed to its neighbors. If 
a = 1/4 the system is conservative. If one of the neigh- 
bors reaches the threshold the process is repeated until 
all the sites have their values below the threshold. The 
number of slips (or the number of sites involved in the 
process) is defined as the avalanche size. Then again, the 
site closest to the threshold is found and a new avalanche 
is triggered. The avalanche distributions follow power 
laws and by varying the degree of dissipation a, the slope 
of the distribution can be tuned. Although the criticality 
of the system in the nonconservative regime (a < 1 /4) 
has been widely debated 0, 0| , the model results in a 
power law distribution of avalanches (for a ~ 0.2) simi- 
lar to the Gutenberg-Richter law 0] and als o rep roduces 
other characteristics of real earthquakes 0, Il6j . 

In the aim of improving resemblance with the mechan- 
ical model we have introduced two variations in the OFC 
model. First: Thresholds are distributed randomly fol- 
lowing a Gaussian distribution of standard deviation a. 
When a block slips a new threshold is imposed to its site. 
This is closer to the actual block-surface friction problem 
and allows the system to start from a configuration of 
zero force in every block [l7|. Second: Instead of assum- 
ing infinitely accurate tuning, we add a quantum of force 
in each step. In the mechanical model this is equivalent to 
a finite velocity (considering finite time resolution). The 
dynamics become more complex because several sites can 
reach the threshold as the plate moves, so several clusters 
start to grow and eventually they can touch each other 
merging into a single one. The avalanche size is defined 
as the number of slips (or number of sites involved) in 
each cluster. 

In order to study the temporal series of avalanches we 
need to define a unit of time. Let us consider the case 
a = 0, a = (isolated blocks with single threshold). 
When the force that drives the system reach a value F m 
equal to the threshold, every site returns to its original 
value and a trivial periodicity rules the dynamics; this 
period is a natural unit. Therefore our unit of time is 
the real time T m that a mechanical model spends to add 
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FIG. 1: Avalanche size distributions for the number of slips 
in each cluster, with a = 0.2 and a = 0.001. The slope 
is —1.91. Inset: collapse of all the curves under the scaling 
relation P{s,L)L fi = f{sL~ v ), with (3 = 4.2 and v = 2.2. 



a force F m to the system. 

We performed simulations taking a quantum of force 
5F = 10~ 4 and threshold values around 1 for different 
a values: 0, 0.001, and 0.01. Then an isolated block 
would need 10 4 steps to reach the threshold, so it would 
spend a time 6t/T m = 10 -4 in each step. We find that 
the avalanche size distributions follow power laws for dif- 
ferent values of a, with cut-offs sensitive to the size of 
the system. The avalanche time series show quasiperi- 
odic behavior with a period proportional to the degree 
of dissipation (similar behavior can be found in the OFC 
model); periodicity is less robust than criticality in the 
system, and foreshocks and aftershocks are strongly re- 
lated with the observed periodicity. All the simulations 
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presented in this paper took place with open boundary 
conditions, but the same results were obtained in a sys- 
tem with free boundary conditions 

Figure 1 shows the avalanche size distribution for dif- 
ferent sizes of the lattice (more than 3 • 10 8 clusters for 
L=512), where avalanche size is defined as the number 
of slips in each cluster. The system, with a = 0.2 and 
a = 0.001, has been driven adding a force SF = 10~ 4 
in each step. The distributions follow a power law with 
a slope equal to -1.91. The curves present cut-offs sensi- 
tive to the size of the system but they collapse (see inset) 
when the finite size scaling relation P(s, L)L@ = f(sL~ y ) 
is applied with (3 = 4.2 and v — 2.2. This v value 
larger than the system dimension indicates that for larger 
avalanches, sites slip more than once. Simulations for 
other a values with L=512 displayed that the absolute 
value of the slope of the distributions increases as the 
dissipation decreases, in contrast to Also the power 
law behavior of the avalanche distributions for small a 
values is more robust in this finite velocity model than in 
the original OFC one due to the fact that larger clusters 
have less probability in the OFC model. The avalanche 
size distributions (at least for a > 0.1) do not suffer con- 
siderable variations when a moves from to 0.01. For 
a = 0.1 the system is not critical anymore, but this is 
not a realistic value for the fluctuations of the friction 
force associated to an interface between a block and a 
flat surface. 
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FIG. 2: Auto-correlation function for the avalanche time 
series (eq:l) for a system with a = 0.001 and L=128. Notice 
in gray the same function for a system with L=128, a = 0.125, 
but with a = 0. 



FIG. 3: above: position of the peaks for the correlation func- 
tion of the avalanche time series (eq:l) for different a. They 
follow the equation T — Tm(l — 4a). below: height of the 
peaks. 



3 



C/5 
CD 

o 

03 
CD 
> 

cc 



CD 



0.9 



0.8 



0.7- 





a = 0.25 





-0.2 0.0 
t 



a = 0.20 



06 ^LAjJukiAM m^UkLiX^liim \l.u 

. ifPPfWr ifff^f^f^f ^ppppfiF *m 



-0.2 0.0 0.2 

t 



FIG. 4: Number of avalanches in the vicinity of a large 
avalanche, normalized to the total number of large avalanches, 
in a system with L=128, a = and a = 0.2. Inset: in a sys- 
tem with L=128, cr = and a = 0.25 (conservative case). 



The analysis of the auto-correlation function 
E«r)- S (T + t))-< S (r) > 2 



C(t) 



£( S (r)- < s(t) >y 
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where s(t) corresponds to the avalanche time series (Fig- 
ures 2 and 3) displays a strong correlation between 
avalanches. Some of the peaks for a = 0.001 are shown 
in Figure 2. The position and height of all the peaks 
for different a appear in Figure 3. The position of the 
peaks indicates that for every analyzed a the system has 
a quasiperiodic behavior with a period proportional to 
the degree of dissipation: 



T m (l - 4a) 



(2) 



In general the height of the peaks decreases as the dis- 
sipation decreases; and for a = 0.225 the only system 
where a sign of periodicity can be seen is for a = 0.001. 
Nevertheless, the peak is just a little higher than the 
background noise (see Figure 2). For the conservative 
case, as in the "sandpile" model, avalanches are uncorre- 
cted. Periodicity in the system is not as robust as criti- 
cality when a varies: peaks are almost delta functions for 
(j = (see curve in gray in Figure 2) but their width in- 
creases dramatically when a increases. The height of the 
peaks shows a monotonous variation with a for a = 0, 
but for larger a there are local variations not well un- 
derstood yet. Although the curves from smaller a val- 
ues cross each other, the peaks for larger a are extremely 
wide, noisy, and have a very small height, indicating that 
periodicity vanishes when noise is added to the system. 

The behavior of the earthquake time series in the vicin- 
ity of a large event have been studied for more than a cen- 
tury, and Omori's law is one of the most well known issues 
in seismology: the rate of aftershocks decays following a 
power law with slope -1 after a main shock. We have char- 
acterized the avalanches in our model into small, medium 




FIG. 5: Foreshocks (left) and aftershocks (right) for the 
avalanches displayed in Figure 4. They have been plotted in 
log-log scale for comparative purposes. The slopes are -0.03 
and -0.05 for foreshocks and aftershocks respectively. 



or large in the following way: In the avalanche size dis- 
tribution, in a log-log plot, we take the linear zone; in 
the case of L=128 it spans from 1 to 6 • 10 3 (see Figure 
1). Then this interval is divided in three zones logarith- 
mically cquispaccd. As a result, avalanches smaller than 
18 are considered small, those lying between 19 and 330 
are medium, and those greater than 330 are large. The 
number of avalanches larger than zero around a large 
one (normalized to the total number of large avalanches) 
for a system with L=128 is displayed in Figure 4. For 
a = 0.2 we see clear signs of periodicity with a period 
equal to 0.2, and the decays in the avalanche rate around 
a large event are strongly connected with this observed 
periodicity. For the conservative case (a = 0.25) the 
large avalanches are not preceded nor succeeded by any 
event in the avalanche time series. If we plot both the be- 
fore and after decays for a = 0.2 in a log-log graph (just 
for comparative purposes) and we fit lines (Figure 5) the 
slopes are -0.03 and -0.05 respectively for foreshocks and 
aftershocks. If instead of the analysis of the total num- 
ber of avalanches around a large shock we do not consider 
the small ones, and the same analysis is made, then both 
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FIG. 6: Power spectrum for the avalanche time series for 
the OFC model (left) and for our model (right). Both with 
L=128, a = 0.15 and a = 0. 
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slopes are around -0.18. These small values are in agree- 
ment with the previous studies in the OFC model |15|. 
but due to the criterion chosen to define a large event, 
they are sensitive to the system size. 

In order to construct an avalanche time series for the 
OFC model we need to define a scale St to measure the 
time. In the finite velocity model it is natural to choose 
5t/T m — SF, we will use the same value for the OFC 
one (the qualitative result is independent of St). If more 
than one event took place in the same St, the value of the 
avalanche size in this interval is equal to the sum of the 
avalanche size in all those events. The power spectrum 
for the avalanche time series for the OFC model with 
L=128 and a = 0.15 appears in Figure 6. It displays a 
sharp peak at a frequency equal to 2.5 corresponding to 
a period equal to 0.4. The power spectrum for a finite 
velocity model with the same L, a, and a = shows no 
qualitative differences with the OFC one, but the inten- 
sity of the peaks are larger. This shows that periodicity 
is not confined to this particular finite velocity model but 
it is a general characteristic of the OFC one. Neverthe- 
less, in this case, as v — > 0, the unit of time T m — > oo. 
The same periodicity was found in the analysis of fore- 
shocks and aftershocks. The correlation function method 
did not bring satisfactory results due to poor statistics for 
small values of a (the peaks appear for the respective pe- 
riods, but the background noise shows large fluctuations 
for small a values). 

Due to the relation between dissipation and periodic- 
ity in the model we performed a few simulations where 
the value of a is not constant but randomly distributed 
following a Gaussian centered in a, with standard devi- 
ations cr Q equal to 0.005, 0.01, and 0.02. For a = 0.2 
and the a a values being 0.005 and 0.01, the height of 
the peak in the correlation function decreases more than 
90% and the period widens up to 5%. The avalanche size 
distributions do not suffer considerable variations. Sim- 
ilar results were obtained for a — 0.15. For a a equal to 
0.02 criticality disappears for both values of a. This cor- 
roborates that periodicity is more fragile than criticality 
when noise is added to the system. 

Quasiperiodic signals in earthquake time series has 
been used in an attempt to predict the next main shock, 
but generally with unsuccessful results [21 |a| • This sit- 
uation, in combination with poor statistics and the lack 
of a theory that explains periodicity have created doubts 
about the real existence of those series of quasiperiodic 
events [lflj . Gao et al f2Q] found an annual periodic- 
ity following the 1992 Landers earthquake in California, 
but they suggest seasonal differences in water extraction 
rates, rainfall and barometric pressure as the cause of it. 
Considering the periodicity found in this simple mapping 
of the block-spring model into a nonconservative cellular 
automata, we can speculate that the earthquake's natural 
behavior is a quasiperiodic state and that the variations 
or absence of periodicity is due to changes in the dissipa- 
tive regime and/or in the relative velocity of the plates 



and/or in the amount of energy that can be stored in a 
given zone between two tectonic plates (related to our 
threshold that rules the unit of time) . 

In conclusion, we have introduced two variations in the 
OFC model improving resemblance with the spring-block 
one, and bridging the gap between the model and possible 
experiments. We found quasiperiodic behavior in the sys- 
tem with a period proportional to the degree of dissipa- 
tion. Foreshocks and aftershocks are strongly connected 
with the observed periodicity, and for small variations in 
the thresholds or in the degree of dissipation, periodic- 
ity tends to vanish, while the system remains showing 
avalanche size distributions that follow power laws. 
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